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By H. Schlichting 


SUMMARY 


In a previous report by W. Tollmien, the stability of laminar flow 
past a flat plate was investigated, by the method, of small vibrations, 
and. the wave le ng th X => 2«/a, the phase velocity c r , and. the Reynolds 

number R of the neutral disturbances established. In connection with 
this, the present report deals with the average disturbance veloc- 
ities \| u* ^ and and the correlation coefficient 

as function of the wall distance y for two special neutral disturbances 
(one at the lower and one at the upper branch of the equilibrium curve 
in the aR plane) . The rami mum value of the last two quantities lies 
in the vicinity of the critical layer where the velocity of the basic 
flow apd the phas e velocity of the disturbance motion are equal. The 
energy balance of the disturbance motion Is investigated. The transfer 
of energy from pr imar y to secondary motion occurs chiefly in the neigh- 
borhood of the critical layer, while the dissipation is almost completely 
confined to a small layer next to the wall. The energy conversion in 
the two explored disturbances is as follows: In one oscillation period, 

half of the total kinetic energy of the disturbance motion on the lower 
branch of the equilibrium curve Is destroyed by dissipation and replaced 
by the energy transferred from the primary to the secondary motion. For 
the disturbance on the upper branch of the equilibrium curve, about a 
fourth of the kinetic energy of the disturbance motion is dissipated and 
replaced in one oscillation period. The requirement that the total 
energy balance for the neutral disturbances be equal to zero Is fulfilled 
with close approximation and affords a welcome check on the previous 
solution of the characteristic value problem. 

*"Amplitudenverteilung und Energiebilanz der kleinen Storungen bei 
der Plattens tromung . " Nachrichten von der Gesellschaft der Vissenschaf ten 
zu Gottingen, Neue Folge, Band 1, No. 1935* PP • 7®» 
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1. Introduction 

The numerous efforts, within the last decade, to solve the problem, 
of turbulence (reference l) have at least produced some satisfactory 
results for a certain class of boundary— layer profiles when their 
stability was investigated by the method of sma.l 1 vibrations, with due 
consideration to the friction of the fluid and the profile curvature 
(reference 2 ). Referring to Tollmien (reference 3), who treated the 
example of laminar flow past a flat plate, the writer investigated 
several other cases: the Couette flow (reference 4), the amplification 

of the unstable disturbances in plate flow (reference 5 ) , and the 
stabilizing effect of a stratification by centrifugal forces (reference 6 ), 
and temperature gradients (reference 7) • Every one of the investigations 
was restricted to the solution of the corresponding characteristic-value 
problems, without calculating the characteristic function itself. In 
that manner, the wave lengths of the unstable, hence "dangerous”, 
disturbances were identified as function of the Reynolds number. In 
most cases, only the disturbances situated right at the boundary, 
between amplification and damping, were determined. For these, just 
as much energy is transferred from the primary to the secondary flow, 
as sec ondary-mot i on energy is dissipated by the friction so that the 
total energy balance is zero. 

All the stability studies made up to now were, for reasons of 
mathematical simplicity, based upon an assumedly plane fundamental flow, 
which depends only on the coordinate transverse to the direction of the 
flow, and a plane superposed disturbance motion which propagates in form 
of a wave motion in the primary— flow direction. While there 1 b no 
objection to the limitation to the plane fundamental flow, since it is 
frequently realized experimentally, objection may be raised to the plane 
disturbance motion because the disturbances accidentally produced in 
practice are almost always three-dimensional. Accordingly, it might 
appear as if the limitation to two-dimensional disturbances was all 
too special. However, H. B. Squire (reference 8 ) recently demonstrated 
on the Couette flow — this theorem is equally applicable to boundary- 
layer profiles — that precisely the specific case of the two-dimensional 
disturbance motion is particularly suitable for the stability study in 
the following sense: According to Squire, a two-dimensional flow, which 

is unstable against three-dimensional disturbances at a certain Reynolds 
number, is unstable against two-dimensional disturbances even at a lower 
Reynolds number. The two-dimensional disturbances are therefore "more 
dangerous" for a flow than the three-dimensional. The critical Reynolds 
number, which is defined as lowest stability limit, is thus obtained 
precisely from the two-dimensional, not the three-dimensional, disturbances 

To gain a deeper insight into the mechanism of the turbulence 
phenomena from small unstable disturbances, a more detailed knowledge 
of the properties of these small disturbances is necessary. The present 
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report, therefore, deals first, with the distribution of the amplitude 
of the disturbance over the flow section, that is, the calculation of 
the characteristic functions and second, with the study of the energy 
distribution and energy balance of the disturbance motion. The investi- 
gations are based upon the disturbances of the 1 aminar flow past a flat 
plate which are situated exactly at the boundary between amplification 
and damping (neutral oscillations). 


Chapter I 

AMPEETUDE DISTRIBUTION 

2. Discussion of the differential equation of disturbance. 

Let U(y) be the velocity distribution of the fundamental flow 
(fig. l) and ijr the flow function of the superposed disturbance motion, 
which is assumed as a wave motion moving in the x direction (direction 
of primary flow) , whose amplitude cp is solely dependent on y, hence 

ilr(x,y,t) = <p(y)e i ^ a,£ “^ t ^ = q>(y) e ia ^ x-ct ^ 

a is real and A. = 2 it /a is the wave length of the disturbance; 

P = 3 r + iPj^ and c = P/a, are, in general, complex; T = 211/3^ is 

the period of oscillation; P 1 indicates the amplification or the damping, 
depending upon whether positive or negative; c r = P r / a is the phase 

velocity of the disturbance . For the disturbance amplitude cp, after 
introduction of dime nsionless variables from the Navier-Stokes differential 
equations, it results in a linear— differential equation of the fourth 
order, the differential equation of the disturbance 

(U - c)(<p" - arcp) - U"cp = - j~-(cp"" ~ 2a?y" + cAp) (l) 

(R = UjjjS/v = Reynolds number, U m = constant velocity outside of the 

boundary layer, & = characteristic length of the boundary— layer 
profile = boundary— layer thickness, v = kinematic viscosity.) The 
general solution cp of the disturbance equation is built up from four 
particular solutions cp^_, . . . , cp^ 

q> = Cj^cpj + CgCPg + ^3^3 + 


( 2 ) 



k- 
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The boundary conditions qp = <p* =0 for y = 0 and y = » give 
(as explained in the report cited in reference ( 5 )) 




= 0 


(3a) 


and for Cl, c 2 , C 3 the system of equations 

C1CP10 + C 2 cp 20 + C 3 cp 3Q = 0 

^1^10* + C 2 Cp 20* + C 3^30* = 0 
C l $ la + C 2 4 ’2a = 0 


w 


. (^ya “ ^va* + “‘Pya' v | 

The subscript 0 indicates the -values at the wall y = 0, subscript a 
the values in the connecting point y = a to the region of constant 
velocity. From (h) the equation of the characteristic value problem 
follows as 


’io 

^20 

^30 



) * 

10 

-9 

ro 

0 

V 

= 0 

( 5 ) 

'la 

®2a 

0 

. 



This equation is discussed in the earlier reports for several cases. 

It contains, aside from the constants of the basic profile, the 
parameters a, E, c , and c^. The complex equation ( 5 ) is equivalent 

to two real equations, and, if limited to the case of neutral disturbances 
(cjl = 0), these two equations give, after elimination of. c r , one 


equation between a and E. This is the equation of the neutral curve 
in the ocB plane, which separates the unstable from the stable disturbance 
attitudes, and was originally computed by Tollmien for the plate flow. 

They are assumed to be known for the subsequent study (fig. 2).^- 


-kfhe writer computed the neutral curve for the plate flow cited 
under reference 5 again and found Borne differences with respect to 
Tollmien' s report. The newly obtained values are used in this study. 
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In the following, the amplitude and. energy distribution is computed 
for two neutral oscillations, one of which lies on the lower, the other 
on the upper branch of the equilibrium curve (fig. 2) „ The parameters 
of these two neutral oscillations, obtained from the earlier calculation, 
are indicated in table 1. 

For the calculation of the integration constants Cj_, C 2 , 
from (10, we put 


C 


1 


= 1 


(3b) 


because the amplitude of disturbance remains indeterminate up to a 
constant factor, the intensity factor of the disturbance motion. Thus, 
for the other two constants 


C 2 ~ 


la 

2a 


( 6 ) 


1 


1 

’ia, , i 

^30 

? _<p 

®2a ?20 * 10 J 

1 ©TT 1- ‘ 

*30 

*2a ?20 10 

*■» 


The particular solutions cp 1 (y) and cp 2 (y) are readily obtainable 

by expansion in series from the so-called frictionless differential 
equation of disturbance of the second order, which follows from the 
general equation (l) by omission of the terms on the right-hand side 
afflicted with the small factor l/aR. The point U = c r , where phase 

velocity of disturbance motion and basic flow velocity agree, and which 
is termed the critical point y = y^, is a singu l a r point of the 

frictionless differential equation of disturbance, which plays a 
prominent part in the investigations. 
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Putting 

y -yk 

y^ 


' y i ; = V ^ - “i ; “ (a - V ' “2 


vhere subscript k denotes the values at the critical points, 
frictionless solutions for linear^ velocity distribution read 


qpi = cc. 1 sinh(aQ_y-]_) i cp = coshCo^y-L) 


and for parabolic velocity distribution 


<P 


a-^ - Va + V2 2 + a 3 y 2 3 + •" 


u H 

cp 2 = b Q + b 1 y 2 + b 2 y 2 2 + ••• + ^7- cp 1 log y 2 for y 2 > 0 

k 


V 


9 2 = h Q + h x y 2 + bgjr^ +••• + — ^Cl°g|y 2 [ ~ for 7 2 


^The Blasius profile of the plate flow is approximated by 
and a quadratic function (fig. l) , namely 

0 ^ y /5 ^ O.I75: U/l^ = 1.68 y /5 

0.175 ^ y/s = 1.015: U/U* = 1 - (1.015 - y/ 5 ) 2 


y/S > 1 . 015 : U /CT m = 1 

For the connection between 6 and the displacement thickness 
in figure 2, 6* = 0.34lS is applicable. 


( 7 ) 

the 

( 8 ) 


( 9 ) 

< 0 


a linear 


( 7 a) 


5 * used 
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According to earlier data, the coefficients are given by the equations 


= -1: 


a, 




a l = I! ^ = - p a 3 = -£-S a 4 = “ ^q“~» 


a a - 0.0013 c^ 2 + 0.0083 a^j ag = 0.0024 a, 2 2 + 0.0006 

„ 2 

bQ = 1; b-j_ = 0; b 2 =» — 1 + -g— (10) 

b 3 = 0.125 + O.O56 a^ 2 ; b^ = 0.021 - 0.14-1 a 2 2 + 0.042 a^ 
b^ = 0.005 + 0.005 a- 2 2 + 0.004 

bg = 0.0015 + 0.0012 c^ 2 - O.OO38 a 2 ^ + 0.0014 a 2 6 


The particular solution cp_^, with its derivatives, is regular 
throughout the entire range of flow (—1 y^_ = 0, 0 ^ y 2 , and can 

be numerically computed with these data. But the particular solution cp 
has a singularity, in which cp £ ’ becomes logarithmically infinite in 
the critical layer y = y^. ThS more detailed discussion haB shown that 

the friction at the wall* and in a restricted vicinity of the critical 
layer, must be taken into consideration. The first gives the friction 
solution cp^, the other, the. friction correction for cp 2 « Introducing 

the new variable 




n - (y - y k ) 



1/3 m 7 ~ y k 

€ 


> 


(H) 


gives (only the greatest terms from (l) for cp(n) being taken into 
account) the differential equation 


iq>"" + t ft" - e 3l 


( 12 ) 


from which follows the correction for cp 2 near the critical layer, as 
well as the third friction— affected solution cp^. The calculation of 
these two solutions merits a little closer study. 



8 


MCA TM 1265 


3. The friction solution cp^ 

The friction solution cp^ obtained from the differential 

equation ( 12 ) when the inhomogeneous terns encumbered with the, small 
factor e are omitted; hence, from the differential equation 

iqp^"" + T)cp 3 " = 0 

Unusual in this equation is that, in contrast to the complete disturbance 
equation (l) and to the frictionlese equation of disturbance, the 
dependence of the parameters a, E, and tl* by (ll) enters only as 
scale factor for ' y, and that it is not at all affected by U". As a 
result, <P 0 (tj) can be computed once for all entirely independent from the 

velocity profile. In this instance, Tiet Jens’ report (reference 9) 
constitutes . a valuable support. A fundamental system cp^” = 

and q> " = F 2 " of equation ( 13 ) is given by 




and t) 


1/2 



1,3/2 e-«A 


(14) 


The expansion in series of the Bessel functions 



v»r(p + v * 1) 
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(l\l/3 1/6 fk ' 

gives , when the constant factors 1 (— i) /rf j 

and f ^(— i) ^ /rf are omitted 


V (t »> - H - 


,13 


19 


2l3 2 *4*7 1»-I3 1, '*1|-*7*10*13 613^*4*7*10 *13*16*19 


+ — 


+ i 


10 


,16 


vw = 1 - 


113 1 *!*. 3'.3^*4*7*10 5l3 5 *li-'7*10*13*l6 

.6 


— + 


y 


> (15) 


.12 


2’.3 2 *2-5 V.3^*2*5-8*ii 


— + 


+ H 


,15 


1 , .3 1 *2 3'.3 3 *2-5-8 5l3^*2*5*8*u*l4 


_ + 


J 


Owing to the boundary condition cp^" = 0 for T = + 00 only the solution 
aggregate 0 -jF^" + $2^2' approaching zero for great positive real t} 

(^1^ = integration constants) comes into question. For great this 

can he represented by the Hankel function of the .second hind with the 

subscript j and the argument e — hence by 


1/2 w ( 2 ) 


9 3 " - “ ,1 S x/3 


!„ 3/2 e -«A 
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Herewith, the looked-for solution of (13) - the constant factor being 
put as = ^o' ^ or sa ^- e simplicity 3 — is given as 

expansion in power series near tj = 0 

rr " ( p iV + W) a&) 

30 

and as asymptotic expansion for great tj 


§7 - *V (*V ■ !) 

Integrating between the limits tj and tj gives 


9 

9 


3_ 

30 


f 



Fl’drj 



P 2 n dTi + 1 


(16b) 


— = h f f S’l'dtldri + P 2 f P F 2 "dTidTi + tj - TJ - D (16c) 

y 30 U\lJ% U\dVo 

where B and D are additional integration constants dependent 
on t) o ; tj 0 is the tj coordinate at the wall, y = 0; hence by (ll) 


/ / s !/3 

% = "^k/ € = (Ha) 

The boundary condition q =1 at y = 0, that is, tj = t| q , is 
fulfilled by adding the term 1 in equation (l6b). 


3 


-Cq9 


3 Y 30 


is the gliding Bpeed of the frictionless solution at the 


wall. 
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Tietjens computed, the Integration constants P^, P 2 , B from 
the fact that at a point 1} = t^, which lies in the range of validity 

of both expansions, asymptotic expansion and power series expansion of 
the solution of (13) up to and including the third differential quotient 
must agree. The resultant equation system, set up and solved by Tietjens, 
reads 


P-lF-l’" + PgFg'" = BF ^” 1 

(a) 

3^2" + P 2 f 2 " = bf 4 " 

(D) 

ab n b 

(c) 

P-L / ^"dTi + P 2 / F 2 "dT) + 1 = BF^’ 


Bo U% 


P 1 P P I 1 2" dTjdTi + P 2 P P F 2 ,, drjdTj + T) - tj o — D = BFjj. (d) 


The Integration constants Pj, P 2 , and D can also be computed 
by a simpler method than TIet Jens', with the aid of the transition 
formula from Hahkel’s to Bessel’s functions, which reads 1 *' 


S l/ 3 (2)(z) ‘ " 15^73 { e ” l/3 J l/ 3 (z) - J -l/ 3 (z> } 

This obviates the Joining of the two expansions, makes (17) superfluous, 
and (17b) gives the exact values immediately: 

h = -V«/ 6 32/3 ill) sni i. e«A 3 l/ 3 r fi) sin | (18) 
Pi „/p\ V 3 / 3 


rf£ 


= -1.190 + 0.6871 


= 0.789(1 + i) 


^Tollmien, who pointed this out to the writer, had this representation 
as far back as 1929, but summarily took over the data by Tietjens for the 
sake of simplicity. 
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The factors B and. D as function of tj can then also, he indicated. 
1 'o 

explicitly, namely 


Bt - < 1 - e 


■" n > m 4) -t -’fX 


0 V dT i 


) _1 

>■ 




y 


D = + P 1 


i° n tj 


’lo^o 


Po p° P 7 ! 

I ’l" dT ) dT l + g“ / / F 2” dT l dT l 


> 


(19) 


- Bcp 3 (0) 


with 


Bcp^(O) = ipj. 


Table 2 contains the results of a new calculation of Tietjens' value 
carried out by these formulas. The differences from Tietjens' figures are 

[)f[ /IT] 

insignificant. The values of Fj", F 2 ", / F^dr] = Fj_' , / F^dy) = , 

(ft pr\ 

Jq. F i' dT l " Fj, / Fg 1 d.T| = F 2 as function of tj are indicated separately 

as real and imaginary part in table 3. Since, according to ( 15 ), these 
quantities are either symmetrical or antisymmetric functions of t], this 
table can be continued i mme diately according to the negative values of tj. 

For the two neutral oscillations, whose amplitude nnd energy distri- 
bution is to be computed, it is 

t] 0 = -2.63 and t] 0 = -4.05 


The corresponding values of c, according to (lla), are given in table I. 
along with the integration constants B-^ B 2 , D obtained by interpolation 

for these tj q values. 

This takes care of all the data necessary for computing the friction 
solution cp 3 with its first and second derivative as function of n by 

the equations (l 6 a, b, and c) . Table 4 gives the thus obtained values 
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* 3 " V ^3 

of ; ; — as function of ti. The connection between q and y 

V ^30 ^30 

is, according to equation (ll) , 

y = 7 k + eTj 

In all cases, the friction solution cp^ from the wall toward the inside 

of the flow is very quickly damped out; but it still extends a little 
beyond the critical layer for both oscillations. 


h. The friction correction of cp 2 in the intermediate layer 

The second frictionless solution cp Q behaves singularly at the critical 

■g «' 

layer y = y fe , namely through equation (9) as - y k )log(y ~ y fe ) , 

so that cp 2 

V 

as -. 

V fy - y k ) 

From the differential equation (12) , in which only the greatest 
friction terms are taken into account, follows a solution cp 2 modified 

by the friction, which joins the frictionless solution at some distance 
from the critical point. For this purpose cp 2 is expanded in powers 

of the previously introduced small quantity e = (dRU^. 1 ) — ^/3 

9 2 = ^20 + 6 ^21 + ’ • ’ ( 20 ) 

Tgo being chosen equal to unity. From (12) follows the inhomogeneous 
differential equation 


V r 1 

* behaves as |l + log(y — y k )j. and 


cp 2 ” 


behaves 


V" - ^21 



21 


( 21 ) 


for 9, 


with reference to t). 



Ik 
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On account of the very small value of e, tj can assume great values 
even at small y — y^ values. An attempt is made to find such a 

solution q>2l' "^is equation, which for large T), hut small (y — y-^) 
joins on to the fhictionless solution 




1 



y -yjt 


For large 


tj there shall he: 




( 22 ) ) 


The corresponding homogeneous equation appeared earlier in the 
calculation of cp^ (equation (13)). It has the fundamental system F-^" (t]) 

and F 2 "(t]) (equation ( 15 ))- A particular solution of (21) is 


<P 21 "(t,) 


-i 






which can he verified easily by substitution, and the general solution 
of (21) is • 


V"(n> - 




F 2 "dn + c^F-l" 



(23) 


The integration constants c^ and Cg, which can he complex, are 
evaluated from the boundary conditions. The quantity cp P -| " is complex 
and shall join the real frictionless value (equation (22)) for large 
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values of tj. By decomposition in real and imaginary parts, the four 
equations defining the integration constants read 



II 

21 * 


G"(ti) = -^"^li* + * 11 * 2 ?' + F lr" F 2i’ 


+ °lr F lr C li F li" + c 2 r F 2 r U “ c 2 i F 2 i" " \ 


V 


5^211" = S ” ( ^ * " F 2i” F li' + ^Ir’ + ^i^l' - (24) 

k 


+ c 


TTt H 

lr^li 


+ c..,F, " 
li lr 


+ c 2r F 21 


+ ° 2 i F 2 r 


= 0 


for tj = i ^ 1 

From (21) , with the boundary condition (22) , it follows that is 

an anti symmetrical,, and cp pn , " a symmetrical function of tj. Moreover, 
since F^', F^’, F^.", F^" are symmetrical and 'F 2s ' , F 2 i"^ F ir" 

anti symmetrical functions of tj, the following must be true 

c li = c 2 r = 0 ( 2 5a) 

for reasons of symmetry. The other two constants c^, c pi are obtained 
by solving the above equation system for tj = tj^. For the present 
calculation, tj^ = 4 was chosen. The series for the Bessel functions 

are still fairly convergent for tj = 4; but since differences of very 
large numbers occur, the separate terms in (24) must be computed to 
five digits (table 3 ) . For c^ and Cg^ 


Op*. = I. 2852 ; c 2i = O .9373 


(25b) 


V V , 

so that — 7 ? pi r ." = G" (tj) and — 7^21i" = can be calculated. The 

k U k 

values are given in table 5 . 
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dc ?2 

The Tallies of cpg ajld ^e dn ^ e:mLedda 't'® layer are obtained 

immediately by quadratures, namely 


and 


f*2r 

ay 


V 


- j^(l + log e) + G* 


(2oa) 


d( P21 

dy 


V*,, > 3l. 

XV s (t > } “ n • 



H"(ri)dTi 



(26b) 


A check on this numerical calculation 1 b given hy the fact that 
d9>2i 

for — at transition from large positive to large negative T) the 

transition substitution for cp 2 deduced hy Tollmien (reference 3) must 

result again (compare equation { 9 )), which he obtained hy discussion of 
the asymptotic representation of the Hankel functions. Tollmien* s 
transition substitution gives 

Ki") . fe gi 

\ dy /y=+°° V dy 



the present numerical calculation gives 


Alep 


2i 


Alep, 


2i 


U " n+00 

+ — / H"(n)drj 


dy /y=+“ \ dy /y=-°° U k * u-00 


and the graphical integration gives 


d-n H" (tj) dTj = 3.1A. (27) 

that is, complete agreement within the Bcope of mathematical accuracy. 
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The intermediate layer near the critical point, "which "by the 
present calculation reaches from about q = — 4 to T) = +4, is already 
so wide at the first neutral oscillation that it reaches up to the "wall 
(wall q = - 2 . 63 ); at the second oscillation with q Q = -4.05, the 

boundary of the intermediate layer is reached exactly at the wall. With 
this, all data needed for the numerical calculation of the solution cp 2 , 

corrected by the friction with cp 2 ‘ and cp " , are available. 


5 . The numerical values of the integration constants 


All three particular solutions 



are numerically known. 


To build up the required solution cp from it, the numerical values of 
the integration constants and must be accertained (equation (6)). 

First of all, equation ( 2 ) is rewritten in a more suitable form, namely 


cp = cp x + C 2 cp 2 + C 3 * 


30* 


<P. 


(2a) 


where equations ( 3 a) and ( 3 b) were resorted to and cp^ was replaced by 
the quantity ^ 3/ ^30 * follows Immediately from the numerical 

calculations. Comparison with ( 6 ) gives 


C 


3 


1 





(6a) 


This method of writing has the advantage that the two integration 
constants C 2 and C^’ in (2a) are dependent only on the values of 

the frictionless solutions cp^ and cp 2 , hence are relatively simple 

to compute. 

The values of cp 1& , cp 2a > 9 la ' j <P 2a ' » ® 2a 811(1 the va l ueB 

° f C 2 811(1 S' thua com P utecl ^ (6) and (6a) for both neutral 

oscillations are given in table 1. Table 6 and figures 3 and k give 
the valueB of cp^, cp 1 , cp^ 1 , cp^ computed with it, hence the desired 

amplitude distribution 'as function of y/6. 
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Outside of the "boundary layer, at y/5 > 1.015 the simple formula 
cp^ = 0*6"°^ cp r ' = -aC*e“ a 5 r 


cp ± = = 0 


( 29 ) 


is valid. for the amplitude distribution. The constant C* is so chOBen 
that the value of <p 1 Joins the already found value in y/5 = 1.015. 
(Table 1.) 


6. The average fluctuation velocities and the correction factor- 

(compare reference 10) 

Changing to the real method of writing 


u' = — = K Jcd ' 


Sy I" 

r « = _ ^ = Ka fi 


cp r ' cos(ax — P r t) — cp 1 ' sin(ax — 3 r t)l 0^ 


(30) 


r t 

^ ~ Ea ^ r sin ^ a ‘ x; ~ + ^i cos(ax — P r t) j- 


U, 


m 


K is a freely available intensity factor. According to figures 3 and 4, 
the one phase (cp r or cp^' ) predominates in both neutral oscillations. 

Tne amplitude distribution of u’ and v' can be represented 
most appropriately by forming, in analogy with the turbulent fluctuation 

\ u' ^ ^ 

velocity, the dimensionless quantities •*— — and 


XL 


XL 


-, where the dash 


m m 

denotes the time average value formation over a period T at a fixed 
point x, y, or in other words 


u 


T 


- u ,2 dt (T = vibration period) 


t=o 
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The result is 


and. 



_E 



+ <Pi 


> 2 . 

f 





(31) 



+ v 


.2 




+ <Pi * 2 + co 2 (cp 2 + 


c Pi 2 ^| 


(32) 


The last quantity gives the mean kinetic energy of the motion 
disturbance. (See eq. 36.) These averages* which are independent of x, 
are represented in figures 5, 6, and 7 and table 7 for both neutral 
vibrations as functions of y/B. The intensity factor itself was so 

chosen that the average value of \u' 2 in the boundary layer is equal 



0.05U. 


mi 


(table l) . 


The maximum amplitude 


for both neutral vibrations lies near the critical layer. The correlation 
factor between u 1 and v* * which is completely independent from the 
intensity of the motion disturbance, can then also be calculated. It is 



k(u* , v‘) 


V<Pi “ Wl* 



+ fPi’ 2 ) (<P r 2 



(33) 


The correlation factor is likewise dependent on y/B only; its 
variation is indicated in figure 8 and table J. It is negative almost 
throughout the entire range of the flow, for both neutral vibrations, 
as is to be expected, since, owing to the positive dU/dy, positive u' 
is usually- coupled with negative v T and negative u* with positive v* . 
The maximum value of k is -O.I7 and -0.19, respectively. It is inter- 
esting to compare the theoretically established correlation coefficient 
with Townend 1 s data in a developed turbulent flow (reference 11). The k 
values of -0.16 to -0.l8, obtained for the flow in a channel of square 
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cross section at various distances from the axis., are of the same order 
of magnitude as those obtained hy the present calculation for the 
incipient turbulence. 


Chapter II 
ENERGY - DISTRIBUTION 

7 . The kinetic energy of the disturbance motion. 


Having established the amplitude distribution for the two neutral 
vibrations, the energy of the disturbance motion can be computed. The 
total kinetic energy of the disturbance motion in a layer of unit height 
which, in x direction, extends over a wave length X an d in y direc 
tion from the wall to infinity is 

P>X Coo 

E = / / (u* ^ + v ,2 )dxdy 

2 </ x=0 \j y=0 


- f ^ 2 u m 2 J Q {jPr * 2 + ^l ' 2 + a2 ($r 2 + c Pi 2 ^j d.(y/ 5 ) (3^) 

The energy dE of the secondary motion in a strip of width dy 
and length X is accordingly 


Besides, 


dE 

ay 


= | |o m 2 K 2 |<p r ' 2 + cp i ' 2 + a 2 (cp r 2 



0.533 


5_ dE 



+ v» 


2 



(35) 


(36) 


E q is the basic— flow energy in a layer of unit height, length X, and 
width S (compare equation (32)). Figure 7 shows the dimensionless 
energy distribution by equations (35) and ( 36). The energy is strongly 
concentrated near the critical layer. To obtain the total energy E, 
the integral (34) must be evaluated. Dividing It in two parts with the 

limits 0 % y/S = 1.015 and 1.015 = y/ 6 < 00 , the first portion is 
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obtained, by graphical integration based, on the computed, amplitude 
distribution. The second portion is obtained analytically by (29), 
namely 


f y/ 8=1.015 L 


9 . 




a 2 (q) r 2 + cp, 2 )] d(y/S) = * 0 ^ a5 


The results of the evaluation are given in table 1. 


Wow the energy of the disturbance motion is compared with the basic- 
flow energy E c in the space of -unit height and surface area 1X6. 

It is, by equation (7a) 


E 


o 




U(y)dxdy = 0.533 



(37) 


Hence, for the ratio of energy of the secondary motion to the energy ' 
of the basic flow E/e o the values presented in table 1 are obtained. 


8. The energy balance of the disturbance motion. 

Consider the time variation of the secondary-motion energy of a 
particle that moves with the basic flow, hence 


D 

Dt 




(38) 


For stable disturbances, the total change of energy of the secondary- 
motion is 2 if 1 ff^ 2 + v '^ ^ ^ or ^stable disturbances > 0, 


and for neutral disturbances = 0, the integration extending over the 
entire range of the particular flow. Participating on the variation 
of the secondary-motion energy are: first, the transfer of kinetic 

energy from the primary to the secondary flow, or vice versa; second, 
the pressure variation; and, third, dissipation. For neutral vibrations, 
the total energy balance is not only eq.ua! to zero for the entire space 
in question, but for every point y of the cross section, the energy . 
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increase per Titration period T = 2 Jt/P r is also eq.ua! to zero. This 
is easily confirmed in the following manner: It is 





u 


+ v ,2 }dt + £■ U 
2 / 


1 T 
t=0 


^y-( u’ 2 + v ' 2 ) dt 
02 





+ v* &-) dt = 0 
3 x dx / 


The first term disappears ty reason of the periodicity of u’ and v' . 
The same holds true when the last term for u* and v* is entered 
according to equation (30) . Thus, the energy increase per vibration 
period T is equal to zero at eveiy point x, y for a neutral 
vibration. 


It is interesting to see how the several factors enumerated above 
participate on the energy conversion In a specific case. For both 
specific cases of neutral disturbance the calculation of the energy 
Is carried out for a plane basic flow and a plane disturbance motion 
according to Lorentz (reference 12 ) 


D 

Dt 


r 


1 



1 2 


+ v' 2 ) 


-pu' V* 


dU 


d(u‘p' ) + dfv'p') [ 



where 




+ P- 






P> = coefficient of viscosity 

ay 
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The first term gives the transfer of energy from the primary to the 
secondary flow, the second gives the contribution resulting from the 
pressure variations, and the third and fourth terms, the loss of energy 
by dissipation. After integration of this term with respect to y over 
the total width of the laminar flow from y = 0 to y = “ and with 
respect to x over a wave length X of the disturbance, the second and 
fourth terms disappear, since u* and v‘ disappear for y = 0 
and y = ■» and with respect to x have the period X. Thus, the 
growth of the energy per unit time in a layer of unit height and base 
area 0<y<oo, 0<x<X is: 


DE 

Dt 


-P 




dU . . 


(y=0 





(39) 


The first integral gives the total energy passing from the primary to 
the secondary motion; the second, the total dissipation. The portion 
of the energy due to pressure variation is removed by the integration. 
The two energy portions for the two neutral disturbances are evaluated. 
Through substitution of (30), followed by integration with respect 
to x, we find 


2 EEL = 
X Dt 


— paK?U 2 
JR 



(<P r ,c Pi - c P r <? i ,) 


SSL 

<Jy 


fly - rtJ m 2 






or 


t ST - ^ 1 ♦ e 2> 


xs £ TJ 2 

2 m 


e^ and eg denoting the dimensionless energy integrals 


(40) 


©i = — 


(cpr’cpi -qypi’) 


*(u/u m ) 


d(y/&) = 


0 


&(y/ 5 ) 


0 


e l’ d (y/ & ) (4la) 
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2 = ~^J Q \^> T " -<£y v ) 2 + (cp 1 ” - a 2 ^) 2 ^ d(y/5) = J e 2 *d(y/S) (4lb) 


To find the energy change of the disturbance motion in the vibration 
period T(T = X/c r , o r = phase velocity), this energy change is referred 

to the total kinetic energy E of the secondary motion which is given 
by equation (34) . From (40) follows then the specific energy change of 
the disturbance motion as 


_T_ DE 
Eq Dt 


U m 2* . , 

c r 0.533 Zj 0 ^ e l 


+ e 2 ’ )d(y/&) 


(42) 


where Z = 0.432 and Z = 0.8l0 for the first and second neutral 
vibration, respectively, while U m /c r =2.86 for both neutral vibrations. 


The local energy transfer from primary to secondary motion (l) and 
the local dissipation ( 2 ) for the two neutral vibrations is then 


5 d 


E 3^1,2 = 78 -° e l,2' 0r = 4l ' 6e l, 2 ' 


(43a, b) 


when 


T 22!. = £E 
Dt 

The values of e ± ’ (equations 4la, b) can be obtained (table 8) on 

the basis of the computed amplitude distribution for both neutral 
vibrations. Figure 9 represents the local energy conversion. The 
dissipation in wall proximity is seen to be extremely great, while the 
critical layer is of no particular importance for the dissipation. But 
the energy transfer from the primary to the secondary motion is greatest 
in the neighborhood of the critical layer, while at the wall and farther 
outside it is very small. The curve is similar to that of the correla- 
tion (fig. 8), as anticipated. 

The graphical integration of e-^ and e 2 ' gives the values 
indicated in table 1. The energy balance = 0, or e-^ + e 2 = 0 
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for the neutral vibrations is therefore fulfilled with satisfactory 
approx imat ion^ , and constitutes a very welcome check on the rather 
complicated solution of the characteristic value problem. 

The total energy transferred in vibration period T from the 
primary to the secondary motion is 

(ae) 1 /e = 78 . 0 ^ or 41.66^ 

and the total energy dissipated 

(AE) 2 /E = 78.0e 2 or hl.6e 2 

These figures are also shown in table 1. Thus, at the first neutral 
vibration, about half of the secondary-motion energy is destroyed by 
dissipation during one vibration; at the second neutral vibration, the 
energy conversion is only about half as great. 

At the second neutral vibration, the vibration period is a little 
greater than at the first, that is, as 1 b readily obtainable from the 
data of table 1, is 


T 1 ■ (g) ‘ 10a X ^ jf? T 2 = (g) 2 ' U ' 8 X 10 

1 m 

To illustrate; For a plate flow in water with 

o 1 1 

V = 0.01 cm^sec , TJ_ = 20 cm sec 


4 V 


U 


m 


^According to the present calculation, the dissipation for both 
vibrations is somewhat greater than the transfer of energy from the 
primary to the secondary flow. This is due to the fact that in the 
stability calculation only the dissipation of the friction solution cp^ 

was taken into account, while the dissipation of the frictionless vibra- 
tion (<£-[_, cp 2 ) was ignored. But, in the energy equation, the dissipation 

of frictionless and frictional vibration was computed and is therefore a 
little greater. Thus, the "neutral vibrations" have, exactly computed, ■ 
still a little damping, and the indifference curve (fig. 2) is, as a 
result, shifted a little toward the inside. 
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the periods of vibration are 

= 2.5O sec; Tg = 3-70 sec 

Thus, vibrations of comparatively great periods are involved. 

Translated by J. Tanier 
National Advisory Committee 
for Aeronautics 
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TABLE 1 


THE PARAMETERS OF THE TWO NEOTRAL VIBRATIONS OF TEE PLATE FLOW 


First neutral Second neutral 



<x5 

(tys/v) 

MV”W> 


^0 

1/e 

Pn 


vibration 


0.466 

2.62 x 10 3 

.163 
= 350 
.209 . 
1.625 
-.494 
-2.63 
12.6 
r .0695 

L+.1021 


r -.1526 
L-.0736 


9i' ? + ct, 2 (cp r 2 + cp ± 2 } J. dy/6 

’O'* 


.07361 

.37^ 

.2001 

.416 


-2.425 

.234 

-2.327 

.101 

1.005 
f -.046 
\+l.563i 
f- 1.00 
-.1571 
.706 
.1454 


Titration 


0*737 , 

6.08 X 10 3 

.258 

.350 

.209 

1.625 

-.494 

-4.05 

19.4 


r -.0470 
l +.02761 


r .0368 

L -.06501 


.435 

.097 

.306 

-2.240 

.417 

-2.014 

.207 

1.011 
r -.114 
1+1.5631 
f —.988 
L-.3251 

1.075 

.1166 


1.015 
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TABLE 1 

THE PARAMETERS OP THE TWO NEUTRAL VIBRATIONS OP THE PLATE PLOW - Concluded 


First neutral vibration 


Second neutral vibration 



0.432K 2 


O.81CK 2 


_E 

E 


0.00913 


O.OllO 


o 

e^ X 10 3 
e 2 X 10 3 

(ae) x 

E 


5-75 

-€.16 

O.kb'j 


6.39 

-7.10 

O.265 


(£E) 

E 


-O.V79 


—0.29^- 
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TABLE 2 


HEWLT CALCULATED VALUES OF ^ AND D 


"^0 

<4 

D 

0 

O.387 + 0.672i 

0.672 - 0.3871 

0.5 

0.341 + 0.3661 

0.770 - 0.3801 

1.0 

0.262 + 0.2131 

0.892 - 0.3501 

1.5 

0.192 + 0,1421 

1.023 - 0.2611 

2.0 

0.132 + 0.1131 

1.202 - 0.1351 

2.'5 

0.0822 + 0.10311 

1.358 + 0.1241 

3.0 

0.0332 + 0.09721 

1.397 + 0.5191 

3-5 

-O.OI65 + O.O782I 

1.139 + 1.0231 

4.0 

-0.0465 + 0.03231 

0.493 + 1.2141 
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TABLE 3 


THE PARTICULAR SOLUTIONS OF cp^ AS FUNCTION OF T) 
[ Compare equation (15) 3 


n 

F lr 


F 

2r 

F 

21 

F, ' 
lr 

F * 
11 

F 1 
2r 

•H 

CM 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.5 

0.021 

0.000 

0.125 

0.000 

0.1250 

0.0005 

0.5000 

0.0026 

1.0 

0.167 

0.003 

0.500 

0.008 

0.4998 

O.OI67 

0.9992 

0.0417 

1.5 

0.563 

0.032 

1.122 

0.063 

1.1186 

0.1264 

1.4865 

0.2065 

2.0 

1.318 

0.177 

1.974 

0.266 

1.9368 

0.5292 

I.8988 

0.6588 

2.5 

2.59 

0.678 

2.97 

0.797 

2.7503 

1.5799 

2.0223 

1.5544 

2.63 

2.87 

0.901 

3.24 

1.022 

2.90 

2.02 

1.95 

1.88 

3.0 

3-97 

1.94 

3.86 

1.91 

2.9208 

3.6995 

1.3354 

2.9269 

3.5 

5.O7 

4.55 

4.02 

3.72 

0.9512 

6.8796 

-I.O85O 

4.2113 

4.0 

4.20 

8.70 

2.35 

5-77 

-5.6379 

9.2445 

-6.1155 

3.3059 

4.05 

3.84 

9.13 

2.00 

5.89 

-6.62 

9.20 

-6.76 

2.90 


antisy. 

ay. 

By. 

antisy. 

By. 

antisy. 

antisy. 

sy. 



0 0 

0.5 O.5OOO 

1.0 0.9980 

1.5 1.4661 

2.0 1.7472 

2.5 1.3100 

2.63 0.941 

3.0 -1.1161 

3.5 -7.6314 

4.0 -19.4902 

4.05 -20.94 

antisy. 


0 

0.0052 

O.O833 

0.4206 

1.3108 

3.0463 

3.64 

5.4735 

6.7189 

0.6805 

-0.806 

sy. 


E 


• 9 999 
.9944 
.9368 
0.6469 
-0.3215 
-0.776 
-2.7429 

-7-2934 

-12.5621 

-12.92 

sy. 


0 

0.0208 

0.1666 

0.5595 

1.2939 
2.3124 
2.58 
3.0209 
1.4618 
—6 . 6661 
-8.10 

antisy . 



O.II87 

-2.28l8 


j. HI 
11 


0 

0.0417 


F 2r’" 


0 

- 0.0010 


0.3331 -0.0333 
1.1165 -0.2525 
2.5541 -1.0533 
4.3761 -3.0889 


Y m 
21 


0 

0.1250 

0.4999 

1.1072 

1.8229 

2.0809 


-8.1612 4.7949 -6.8785 0.1417 
-18.5245 -1.7400 -11.0222 -7.8189 
-27.2703 -26.6992 -7.4580 -26.7961 

ay. antisy. antlsy. By. 
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TABLE 4 


THE FRICTION SOLUTION cpg AS FUNCTION 
OF y/B AND r\ 


First Neutral Vibration 


T1 

f/ 8 

V 

V 

„ tl 

9 3 i 

V 

V 

v 

V 

v- 

V 

*31 

V 

- 2.63 

0 

- 6.40 

7.65 

1 

0 

-o.lll 

-0.0159 

-2 

.050 

- 5.67 

.794 

.665 

.188 

-.069 

-.00 99 

-1 

.130 

- 3.04 

-1.801 

.329 

.106 

-.029 

.0029 

0 

.209 

- 1.92 

-.932 

.136 

-.012 

-.00 9 

.0061 

1 

.288 

-.995 

.113 

.019 

-.040 

-.001 

.0035 

2 

.368 

-.202 

.290 

-.025 

-.018 

0 

.0010 

3 

.447 

.063 

.076 

-.027 

-.003 

0 

.0005 

k 

• 577 

-.O76 

-.012 

-.020 

.012 

0 

.0006 


Second Neutral Vibration 


-4.05 

0 

-17.7 

34 

1 

0 

- 0.0204 

-0.0639 

-3.5 

.029 

-17.6 

8.22 

.460 

.519 

.0001 

-.0554 

-3 

.054 

-9.71 

-3.08 

.106 

.534 

.0070 

-.0415 

-2.63 

.074 

-4.91 

-4.67 

-.030 

.455 

.0076 

-.0322 

-2 

.105 

-.272 

-3.88 

-.103 

.309 

.0051 

-.0202 

-1 

.157 

1.378 

.-2 

-.058 

.166 

.0005 

-.0091 

0 

.209 

.717 

-1.260 

0 

.084 

-.0008 

-.0035 

1 

.261 

-.019 

-.679 

.016 

.035 

-.0002 

-.0012 

2 

.313 

-.194 

-.155 

.007 

.013 

0 

0 

3 

.361 

-.058 

.019 

0 

.012 

0 

0 

4 

.415 

.039 

0 

-.005 

.014 

0 

0 
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TABLE 5 


THE FRICTION CORRECTION q> 3 IN THE 
INTERMEDIATE LAIER AS FUNCTION OF tj 


1\ 

Or" (Tj) 

tf'(n) 

G’(ti) 


0 

0 

1.285 

-0.774 

-1.570 

.5 

.443 

1.165 

-.659 

-.953 

1 

.746 

.857 

-.354 

-.448 

1.5 

.839 

.483 

.051 

-.118 

2 

.767 

.160 

.458 

.035 

2.5 

.589 

-.018 

.798 

.073 

3 

.438 

-.072 

1.056 

.048 

3.5 

.325 

-.055 

1.243 

.013 

4 

.250 

0 

. 

1.386 

0 


&"H) = -a"(n)j h"H) = H"( n ) 

Gr’ (tj) = Gr’ (tj) J H' (-T)) * -Ot - H* (tj) 


\ 
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TABLE 6 


DISTRIBUTION OF AMPLITUDES cp^, Cp ± , <p ' , cp ± * , <p ", qp " 

AS FUNCTION OF 7/ 6 
First Neutral Vibration 


y/ 5 

?r 

<Pi 

V 

V 

cp " 
v r 

’i" 

0 

0 

0 

0 

0 

7.885 

-6.615 

.050 

.Oil 

-.0040 

.411 

-.135 

6.235 

-.003 

.090 

.032 

-.0080 

.620 

-.115 

4.580 

• 770 

.130 

.059 

-.0112 

.782 

-.025 

3.200 

1.717 

.170 

.090 

-.0111 

.896 

.042 

2.355 

1.240 

.209 

.127 

-.0088 

.976 

.069 

I.778 

.434 

.250 

.1 66 

-.0063 

.986 

.070 

-.112 

-.220 

.290 

.203 

-.0041 

.975 

.058 

-.699 

-.576 

.370 

.276 

-.0010 

.831 . 

.018 

- 1.543 

-.347 

.451 

.335 

0 

.670 

.005 

- 1.778 

-.015 

• 531 

.380 

0 

.510 

-.003 

- 1.641 

-.016 

.612 

.414 

0 

.361 

0 

-1.606 ~ 

0 

.693 

.438 

0 

.237 

0 

- 1.514 

0 

.77^ 

.453 

0 

.118 

0 

- 1.431 

0 

.854 

.458 

0 

.007 

0 

- 1.364 

0 

• 935 

.455 

0 

-.101 

0 

- 1.303 

0 

1.015 

.445 

0 

-.205 

0 

-1.252 

0 
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TABLE 6 


DISTRIBUTION OF AMPLITUDES cp r , cp ±J> cp^’ , cp ± * „ cp^", cp ± " 
AS FUNCTION OF y/8 - Concluded. 


Second Neutral Vibration 


y/ s 

<Pr 

<Pi 

V 

V 

cp ” 

*1" 

0 

0 

0 

0 

0 

28.5I 

- 27.89 

.029 

.010 

-.0038 

.720 

-.337 

20.27 

-2.30 

.054 

.033 

-. 01114 - 

1.101 

-.237 

9-16 

6.35 

.07^ 

.057 

-.OllO 

1.215 

-.113 

1.25 

6.28 

.105 

.096 

-.0118 

1.333 

.052 

.I78 

3.60 

.157 

.168 

-.0086 

1.350 

.129 

.558 

-27 

.209 

.236 

-.0011 

1.306 

.076 

-1.218 

- 1 . 5 l 

.250 

.287 

.0005 

1.203 

.oil 

- 2.967 

- 1.31 

.290 

.333 

.0002 

1.091 

-.016 

-2.961 

-.331 

.370 

.409 

0 

.827 

-.016 

-2.86 

.111 

A 5 T 

.1+69 

0 

.607 

0 

-2.58 

0 

• 531 

.507 

0 

.109 

0 

-2.21 

0 

.612 

• 532 

0 

.2I7 

0 

- 1.93 

0 

.693 

. 5 V 7 

0 

.105 

0 

- 1.72 

0 

• 77 l 

• 550 

0 

-.026 

0 

- 1.56 

0 

.851 

.5H 

0 

-.118 

0 

-l.'ll 

0 

.935 

.527 

0 

-.260 

0 

- 1.32 

0 

I.015 

.508 

0 

-.367 

0 

- 1.23 

0 
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TABLE 7 


THE MEAN FLUCTUATION VELOCITIES 



THE KINETIC 


ENERGY OF THE DISTURBANCE MOTION u’ 2 + v' 2 , AND THE 


CORRELATION COEFFICIENT k AS FUNCTION 



OF 

7/6. [equations (31) 

First Neutral Vi' 

, ( 32 ), ( 33 )] 
brat ion 

y/ 5 


u’ 2 

Uv 

cu 

c 

r 

lo2 «' 2 ♦ T ’ 2 

-k 

HI 

- -- - 1 


n 2 

m 

0 

0 


0 

0 

0 

.050 

.445 

.0546 

.198 

.032 

.090 

.647 

.158 

.420 

.061 

.130 

.804 

.287 

.648 

-155 

.170 

.922 

.434 

.854 

.169 

.^09 

I.005 

.608 

1.019 

.140 

.250 

1.016 

• 795 

1.040 

.109 

.290 

1.003 

• 973 

1.019 

.064 

.370 

.854 

1.322 

.748 

.022 

. 1+51 

.689 

1.605 

.501 

.OO75 

• 531 

• 524 

1.820 

• 309 

0 

.612 

• 371 

1.984 

•177 

0 

• 693 

.244 

2.098 

.104 

0 

.774 

.121 

2.170 

.062 

0 

.854 

.007 

2.193 

.048 

0 

• 935 

.103 

2.180 

.059 

0 

1.015 

.211 

2.108 

.090 

0 

1.1 

.203 

2.027 

.082 

0 

1.2 

•193 

1-935 

.075 

0 

1-3 

.185 

1.849 

.068 

0 

1.4 

.176 

1.763 

.062 

0 

1.5 

.168 

1.681 

.057 

0 
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TABLE 7 


THE MEM FLUCTUATING VELOCITIES 



THE KINETIC ENERGY OF THE DISTURBMCE MOTION u ' 2 + v' 2 , 
■ MD THE CORRELATION COEFFICIENT ’ k AS FUNCTION OF y/ 6 . 
EQUATIONS ( 31 ), ( 32 ), ( 33 ) -Concluded 


Second Neutral Vibration 


y/ B 

\L» 2 " 

TO - 

i 02 ^ 

10 2 u’ 2 + v’ 2 

-k 

m 

< 

m 

U 2 
m. 

0 

0 

0 

0 

0 

.029 

.661 

.0655 

•433 

0 

.054 

.937 

.214 

.869 

.121 

.074 

1.038 

.359 

1.068 

.170 

.105 

1.108 

.595 

1.218 

.191 

.157 

1.127 

1.028 

1.262 

.146 

.209 

1.086 

1.445 

1.190 

.076 

.250 

.998 

1.755 

1.019 

.010 

.290 

.906 

2.04 

.854 

-.015 

.370 

.686 

2.50 

.528 

-.019 

.451 

.504 

2.87 

.332 

0 

.531 

.340 

3-10 

.208 

0 

.612 

.205 

3.26 

.146 

0 

.693 

.087 

3-35 

.118 

0 

• 774 

.022 

3.37 

.117 

0 

.854 

.123 

3.33 

.124 

0 

•935 

.216 

3-22 

.149 

0 

1.015 

.305 

3.11 

.187 

0 

1.1 

.292 

2.925 

.171 

0 

1.2 

.272 

2.72 

.146 

0 

1.3 

.252 

2.53 

.126 

0 

1.4 

.234 

2-35 

.109 

0 

1.5 

.218 

2.18 

.093 

0 
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TABLE 8 


THE LOCAL ENERGY CONVERSION l) = TRANSFER FROM 
PRIMARY TO SECONDARY MOTION, 2) = DISSIPATION 
EQUATIONS (41) AND (43) 

First Neutral Vibration 


y/s 

e' X 10 3 

e 2 ' X 10 3 

E dy 1 

1 

0 

0 

86.9 

0 

-6.78 

.050 

.268 

31.8 

.021 

-2.48 

.090 

2.150 

17.6 

.168 

-1.37 

.130 

12.25 

10.8 

.955 

-.84 

.170 

23.25 

5.73 

1.814 

-.45 

.209 

28 

2.64 

2.18 

-.21 

.250 

27.30 

.06 

2.13 

0 

.290 

22.85 

• 72 

I.78 

-.05 

.370 

7*^9 

2.20 

.584 

-.17 

.451 

I.885 

2.80 

.147 

-.22 

• 53i 

-1.102 

2.43 

—.086 

-.19 

.612 

0 

2.36 

0 

-.18 

.693 

0 

2.12 

0 

-.17 

.774 

0 

1.92 

0 

-.15 

.854 

0 

1.75 

0 

-.14 

• 935 

0 

1.61 

0 

-.13 

1.015 

0 

1.50 

0 

-.12 
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TABLE 8 


THE LOCAL ENERGY CONVERSION l) = TRANSFER FROM 
PRIMARY TO SECONDARY MOTION, 2) = DISSIPATION 
[EQUATIONS (la) AND (1+3)} - Concluded. 


Second. Neutral VI oration 


y/s 

e 1 1 X 10 3 

e 2 ’ x 10 3 

13^1 


0 

0 

320 

0 

-13.31 

.029 

1 

83.6 

-.044 

-3.W 

.054 

•7-84 

24.3 

.326 

-1.01 

.074 

18.45 

11.15 

.768 

-.464 

.105 

41.5 

2.61 

1.726 

-.109 

.157 

55-6 

.04 

2.310 

-.002 

.209 

37-5 

.840 

1.560 

-.035 

.250 

5.22 

2.31 

.217 

-.096 

.290 

1.-91 

2.01 

-.331 

-.084 

• 3TO 

-8.44 

1.91 

-.351 

-.079 

.451 

0 

1.80 

0 

-.075 

• 531 

0 

1.38 

0 

-.057 

.612 

0 

1.01 

0 

-.042 

.693 

0 

.819 

0 

-.034 

•774 

0 

.6 95 

0 

-.029 

.854 

0 

.608 

0 

-.025 

.935 

0 

.520 

0 

-.022 

1.015 

0 

.458 

0 

-.019 
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Figure 5.- The mean fluctuating velocity in the x direction 
plotted against the wall distance for both neutral vibrations. 





Figure 6.- The mean fluctuating velocity in the y direction 
plotted against the wall distance for both neutral vibrations. 
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Figure 9.- The local energy conversion of the secondary motion for the 
first and second neutral vibration. I (1), n (1) = energy transfer from 
primary to secondary flow; I (2), II (2) = dissipation. 


NACA-Langley - 4-18-50 - 950 




